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Evaluation  of  the  Magnetic  Induction  Conductivity  Method 
for  Detecting  Frazil  Ice  Deposits 

STEVEN  A.  ARCONE,  BRUCE  E.  BROCKETT,  DANIEL  E.  LAWSON 
AND  EDWARD  F.  CHACHO,  JR. 


INTRODUCTION 

Frazil  ice  is  the  predominant  type  of  ice  found 
in  most  rivers  and  streams.  Information  on  frazil 
ice  formation  is  required  to  understand  the  winter 
hydrology  of  streams  in  cold  regions  as  well  as  to 
plan  navigation  routes,  design  of  bridge  piers,  and 
for  erosion  control,  water  intakes,  and  pipelines. 
Frazil  dams  can  form  under  an  ice  cover,  extend  to 
the  river  bed,  and  impede  or  channel  flow.  As 
there  may  be  no  or  subtle  surface  manifestations 
of  these  frazil  deposits,  time-consuming  drilling  is 
required  to  locate  and  map  them.  Consequently,  a 
remote  sensing  method  was  desired  to  expedite 
mapping.  This  paper  discusses  the  use  of  a  mag¬ 
netic  induction  technique  (MI),  which  is  a  noncon¬ 
tact  method  of  measuring  ground  electrical  con¬ 
ductivity. 

The  initial  evaluation  of  MI  for  locating  frazil  ice 
deposits  was  undertaken  in  1984  as  part  of  a  study 
of  the  effects  of  an  ice  cover  on  the  winter  hydro- 
logic  regime  of  the  Tanana  River  near  Fairbanks, 
Alaska  (Lawson  et  al.  1986a).  Subice  observations 
through  drill  holes  indicated  that  extensive  de¬ 
posits  of  frazil  ice  tended  to  channelize  flow  (Law- 
son  et  al.  1986b).  During  the  1984  studies,  both 
the  MI  and  ground-penetrating  radar  (GPR)  meth¬ 
ods  (see  e.g.  Arcone  [1985]  or  Arcone  and  Delaney 
[in  prep.]  for  description  of  GPR  in  ice  surveying) 
were  tested  to  determine  if  they  could  identify  the 
frazil  deposits  and  channels.  Extensive  data  from 
the  access  holes  were  used  to  evaluate  the  geophys¬ 
ical  data.  The  results  indicated  that  GPR  detection 
of  frazil  ice  deposits  was  questionable  because:  a) 
the  frazil  formations  were  highly  reflective  due  to 
water  content  and  did  not  allow  penetration  (Ar¬ 
cone  and  Delaney,  in  prep.),  b)  the  interface  be¬ 
tween  solid  frazil  deposits  and  the  river  bed  could 
not  be  distinguished,  or  c)  direct  coupling  between 
antennas  offered  too  much  interference.  How¬ 


ever,  good  qualitative  correlations  were  found  be¬ 
tween  the  MI  response  and  the  presence  of  open 
channels  and  frazil  ice  deposits. 

The  objective  of  this  study,  therefore,  was  to 
assess  the  ability  of  the  magnetic  induction  meth¬ 
od  to  map  frazil  ice  deposits  in  an  area  of  heavy 
frazil  ice  accumulation.  Three  parallel  cross  sec¬ 
tions  ranging  from  170  to  210  m  in  length  across 
the  Tanana  River  were  profiled  for  apparent  con¬ 
ductivity  using  MI.  The  results  were  compared 
with  apparent  conductivity  values  generated  by 
theoretical  models  based  upon  extensive  data  from 
drill  holes  to  determine  the  conductivity  of  frazil 
ice  and  river  bottom  sediments  and  the  sensitivity 
of  the  MI  method  to  frazil  ice. 


MAGNETIC  INDUCTION 
CONDUCTIVITY  METHOD 

The  MI  technique  of  conductivity  mapping  util¬ 
izes  the  coupling  between  two  coplanar  antenna 
loops  to  measure  an  apparent  ground  conductiv¬ 
ity:  apparent  in  that  it  is  an  integration  of  any  in¬ 
homogeneities  in  the  ground.  A  ground-based, 
portable  system  has  been  applied  extensively  in 
arctic  studies  (e.g.  Hoekstra  1978,  Arcone  et  al. 
1978,  1979)  where  frozen  ground  or  massive  ice 
formations  can  offer  a  noticeable  conductivity 
contrast  with  the  general  ground  conditions.  Its 
use  for  measuring  sea  ice  thickness  has  been  re¬ 
ported  by  Hoekstra  et  al.  (1979)  and  large-scale 
airborne  systems  have  also  been  tried  (e.g  Fraser 
1978). 

The  MI  method  derives  ground  conductivity 
values  from  the  magnetic  field  coupling  between 
two  loop  antennas  located  at,  or  slightly  above, 
the  earth’s  surface  (Fig.  1).  One  loop,  the  trans¬ 
mitter  antenna,  generatesaprimary,  time-harmon¬ 
ic  magnetic  field,  which  couples  directly  with  the 
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Figure  1.  Schematic  of  the  magnetic  induction  instrument  (Geonics 
EM31)  for  measuring  ground  conductivity.  The  transmitting  loop  T produc¬ 
es  a  primary  magnetic  field  Hp  (which  induces  eddy  currents  Se  within  the  ground. 
Je  then  produces  a  secondary  magnetic  field  Hs,  which  is  received  out  of  phase 
with  H p  at  the  receiver  R.  The  quadrature  phase  component  of  Hs/Hp  is  calibrat¬ 
ed  in  mS/m  of  conductivity.  Also  shown  are  some  of  the  possible  loop  orienta¬ 
tions,  two  of  which — HCP  and  VCP — are  possible  with  the  Geonics  EM  31). 


receiver  loop  but  also  induces  electrical  currents 
(sometimes  referred  to  as  “eddy”  currents)  within 
the  ground.  These  currents  then  generate  a  secon¬ 
dary  magnetic  field,  which  also  couples  with  the 
receiver  loop.  The  secondary  coupling  depends  on 
ground  conductivity.  The  ratio  of  secondary  to 
primary  coupling  is  usually  simply  calibrated 
against  conductivity  so  that  this  ratio  is  the  quanti¬ 
ty  measured. 

Information  about  conductivity  variations  with 
depth  are  achieved  by  varying  either  the  interloop 
spacing  or  the  loop  orientation.  Some  of  the  dif¬ 
ferent  loop  orientations  are  also  shown  in  Figure 
1.  They  are  defined  as  vertical  coplanar  (VCP), 
horizontal  coplanar  (HCP),  and  coaxial  (CA). 
Other  asymmetric  configurations  are  also  possi¬ 
ble. 

The  magnetic  coupling  is  determined  from  the 
mutual  impedance  Z,  defined  as  V/I,  where  V  is 
the  voltage  in  the  transmitter  loop  and  I  is  the  cur¬ 
rent  induced  in  the  receiver  loop.  Z  is  determined 
by  the  linear  combination  of  primary  Vp  and  sec¬ 
ondary  V%  coupling  such  that 


For  this  study,  the  Geonics*  EM31  was  used  in 
the  HCP  mode.  The  EM31  uses  anttnna  loops 
spaced  3.66  m  apart  and  operates  at  39  kHz.t 
These  parameters  and  the  normal  range  of  ground 
conductivity  o  (0. 1  to  10"4  S/m  or  mhos/m)  allow 
the  ratio  Z/Z0  (Keller  and  Frischknecht  1966,  p. 
335)  to  be  approximated  for  HCP  antennas  lying 
on  homogeneous  ground  by  the  formula 


because  the  magnitude  of  yr  is  much  less  than  uni¬ 
ty.  Assuming  negligible  displacement  currents,  the 
quantity  y  is  related  to  o  (the  inverse  of  resistivity 
e)  by  the  formula 

y  =  \iwy*o,  (4) 

so  that  eq  3  becomes 
Z 

-j-  =  1  +  iwpoar/A.  (5) 


V  +  V 
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(i) 


Z  is  normalized  by  the  free  space  impedance  Z0  (in 
effect,  a  measure  of  the  primary  coupling),  which 
is  determined  by  the  formula  (Keller  and  Frisch¬ 
knecht  1966) 


Z, 


ipo^nfA^n^,) 

4irr5 


(2) 


where 

Mo 

i 

n 

r 


frequency  (rad/s) 

4x  x  10*’  (Henry/m) 

number  of  turns 

area  of  the  transmit  ( t )  loop  antenna 
area  of  the  receive  loop  antenna 
loop  separation. 


The  real  part  is  known  as  the  in-phase  component 
and  the  imaginary  part  as  the  quadrature  phase 
component  of  the  normalized  mutual  impedance. 
The  EM3 1  was  used  only  in  the  quadrature  mode 
for  this  study. 

When  the  earth  is  layered,  an  apparent  conduc¬ 
tivity  a,  is  then  defined.  This  quantity  is  equal  to 
the  conductivity  that  would  produce  the  same 
modulus  for  the  mutual  impedance  above  a  ho¬ 
mogeneous  earth.  Interpolation  curves  or  compu¬ 
tational  integrations  of  integral  equations  (Keller 
and  Frischknecht  1966,  Sinha  1976)  must  then  be 
used  to  resolve  the  different  layer  resistivities  and 
thicknesses.  We  used  interpretation  formulas  pub- 
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Figure  2.  Calibration  curve  for  an  EM31  instrument 
held  at  1-m  height  above  an  ice  sheet  of  varying 
thickness  and  a  deep  water  layer. 


lished  by  Geonics  Ltd.  (App.  B)  in  their  operation 
manual  and  based  on  the  work  of  Sinha  (1976). 
Generally,  the  instrument  in  the  HCP  mode  is  sen¬ 
sitive  to  a  depth  of  about  two  times  the  intercoil 
spacing  or  7  m.  Its  lateral  sensitivity  is  about  equal 
to  the  intercoil  spacing.  The  VCP  mode  is  less  sen¬ 
sitive  to  depth.  Arcone  et  al.  (1979)  discuss  this 
more  fully  and  use  a  thin  conductive  layer  at  depth 
to  demonstrate  depth  sensitivity  theoretically.  The 
HCP  mode  was  chosen  because  we  expected  frazil 
deposits  to  6  m  depth  and  also  wanted  to  sense 
bottom  sediments. 

MI  is  most  accurate  for  conductivities  greater 
than  1  mS/m  and  least  accurate  for  conductivities 
less  than  0.05  mS/m  .  This  loss  in  accuracy  results 
from  the  difficulty  in  obtaining  a  site  of  extremely 
low  conductivity  at  which  a  correct  zero  point  (a 
=  0)  calibration  can  be  made.  In  these  studies,  we 
expected  a  wide  range  of  values  due  to  the  general 
four-layer  system  of  ice,  water,  frazil  ice,  and  bot¬ 
tom  sediment,  the  total  depth  of  which  was  within 
the  general  sensitivity  range  of  this  instrument. 
Figure  2  shows  the  calibration  curve  for  solid  ice 
(<x  <  10’’  S/m)  above  Tanana  River  water  (0.01 1 
S/m)  and  gives  values  between  6  and  9  mS/m  for 
ice  of  varying  thickness  over  water  deeper  than 
about  6  m.  Values  less  than  6  mS/m  would  be  ex¬ 
pected  to  have  been  influenced  by  either  frazil  ice 
or  bottom  sediments. 


Figure  3.  Location  of  the  study  site  on  the 
Tanana  River,  near  Fairbanks,  Alaska. 
Cross  sections  of  lines  X6,  X4,  and  X3A  were 
the  3  of  27  transects  that  were  measured  and 
compared  to  MI  data. 


SITE  DESCRIPTION  AND 
SURVEY  METHODS 

The  study  area  is  located  at  the  confluence  of 
the  Chena  and  Tanana  Rivers,  just  south  of  Fair¬ 
banks  International  Airport  (Fig.  3).  The  active 
floodplain  of  the  Tanana  River  in  the  study  area  is 
approximately  500  m  wide,  with  total  water  depth 
ranging  from  1 1  m  in  the  main  channel  along  the 
north  shore  to  about  2  m  over  some  mid-channel 
bars.  MI  apparent  conductivity  and  ground  truth 
data  were  measured  along  three  cross  sections 
(X3A,  X4  and  X6)  as  well  as  along  24  additional 
transects  for  which  ground  truth  was  not  obtained. 
The  EM31  was  hand-held  at  the  operator’s  hip  in 
the  HCP  mode  with  the  instrument  axis  perpen¬ 
dicular  to  the  transects.  Readings  were  taken  every 
2  m  and  typically  30  minutes  were  required  to 
complete  each  profile. 

Hydrographic  cross  sections  of  X3A,  X4,  and 
X6  were  constructed  from  data  measured  in  access 
holes  augered  every  10  m  along  each  cross  section. 
These  data  included  ice  cover  thickness,  frazil  de¬ 
posit  thickness,  and  water  depth.  Snow  depth  was 
also  recorded  to  position  the  instrument  accurate¬ 
ly  above  the  ice  cover.  Data  for  the  frazil  deposits 
describe  only  particle  size  and  sediment  concentra¬ 
tion  and  did  not  record  volumetric  water  content, 
which  was  estimated  to  be  as  high  as  60% .  A  typi- 
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cal  well  log  description  is  given  by  Chacho  et  al. 
(1986).  Water  conductivity  was  measured  at  0.01 1 
S/m  using  a  YSI  conductivity  meter  in  the  field. 


CROSS  SECTION  FIELD  DATA 
AND  MODELING  RESULTS 

The  theory  used  to  develop  the  model  assumes 
homogeneous,  laterally  infinite  layers,  each  char¬ 
acterized  by  one  conductivity  and  the  thickness  re¬ 
corded  in  the  well  log.  In  other  words,  at  each  ob¬ 
servation  point  each  material  (frazil  ice,  solid  ice, 
water,  or  sediment)  is  idealized  as  being  of  infinite 
lateral  extent.  This  is  a  fairly  good  assumption 
considering  the  large  vertical-to-horizontal  scale 
(aspect  ratio  of  ~  20:1)  distortions  used  in  the 
cross  sections  and  in  light  of  the  fact  that  the  in¬ 
strument  is  sensitive  mainly  to  a  lateral  distance  r 
from  the  center  of  the  instrument  (Arcone  et  al. 
1979).  The  modeling  also  accounts  for  the  height 
at  which  the  instrument  was  held.  Conductivity 
values  for  air,  snow  (0  to  35  cm),  and  solid  ice 
were  assumed  10~’  S/m;  any  smaller  value  does 
not  affect  the  results.  Conductivity  values  for  the 
frazil  ice  and  river  bed  were  derived  by  trial  and 
error  for  cross  section  X6  and  then  applied  to  X3A 
and  X4. 


All  field  data  were  gathered  on  1  and  2  March 
1986,  when  air  temperature  never  exceeded  -  2°F. 
This  fact  and  the  diagreement  between  measured 
and  theoretical  values  over  the  deep  water  chan¬ 
nels  on  cross  section  X4  suggest  that  we  may,  at 
times,  have  experienced  operating  difficulties  with 
the  EM31.  Sources  of  error  are  discussed  in  Ap¬ 
pendix  A. 

X6 

The  measured  cross  section  is  compared  with 
measured  (solid  curve)  and  hypothetically  modeled 
(broken  curves)  MI  conductivity  profiles  in  Figure 
4.  The  highest  measured  conductivity  values  oc¬ 
curred  over  the  deepest  water  channels  beneath  the 
ice  cover,  while  the  lowest  values  were  measured 
on  the  north  bank  of  frozen  silt  and  over  the  ice- 
covered  gravel  bar  at  the  southern  end  of  the  cross 
section.  A  large  minimum  occurred  between  50 
and  80  m,  above  the  frazil  deposit.  A  secondary 
minimum  occurred  at  about  125  m,  over  a  small, 
hanging  frazil  deposit. 

A  value  of  1.0x10'’  S/m  was  derived  for  the 
river  bed  by  modeling  the  MI  response  at  points  A 
and  C.  This  value  was  then  applied  at  point  B  to 
derive  a  value  of  6.3  x  10‘4  S/m  for  the  frazil  ice 
deposit.  Ground  conductivity  values  calculated  or 
measured  are  given  in  Figure  4  and  were  used  to 


Figure  4.  Comparison  of  measured  EM 3 1  data 
for  cross  section  X6  (solid  line )  with  a  theoreti¬ 
cal  model  (dashed  line).  Conductivity  values  (S/m) 
for  the  model:  ice  =  10  '  (assumed),  water  =  0.01 1 
(measured),  bottom  sediments  =  10  ’,  and  frazil  ice 
deposits  =  6.3xl0‘.  The  sediment  and  frazil  values 
were  derived  by  trial  and  error  from  data  at  points  A , 
B,  and  C.  The  dotted  curve  is  the  theoretical  response 
if  the  frazil  ice  deposits  were  replaced  by  water.  A 
comparison  of  the  two  broken  curves  demonstrates 
t„e  sensitivity  of  the  Ml  technique  to  the  presence  of 
the  frazil  deposits. 
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produce  the  theoretical  dashed  curve,  which  corre¬ 
lates  well  with  the  measured  curve.  The  dotted 
curve  of  Figure  4  gives  the  theoretical  response 
when  water  is  substituted  for  the  frazil  ice.  It  is  ap¬ 
parent  from  the  comparison  of  the  broken  curves 
here  and  in  the  figures  that  follow  that  the  instru¬ 
ment  is  sensitive  to  the  presence  of  the  frazil  ice 
deposits.  This  sensitivity  reflects  the  low  conduc¬ 
tivity  of  these  deposits,  which  is  surprising  in  view 
of  the  large  amount  of  unfrozen  interstitial  water 
that  was  observed.  Theoretical  mixing  formulas 
(e.g.  Fricke  1924)  suggest  that  ice-water  mixtures 
of  high  volumetric  ice  content  should  be  an  order 
of  magnitude  greater  in  conductivity.  However, 
these  formulations  assume  structures  with  high 
electrical  permeability,  which  may  not  be  true 
here. 

X3A 

The  cross  sections  and  conductivity  profiles  for 
X3  A  are  given  in  Figure  5.  As  in  Figure  4,  the  solid 
curve  on  the  bottom  graph  gives  the  measured 
data,  the  dashed  curve  gives  the  model  data,  and 
the  dotted  curve  gives  model  data  with  frazil  ice 
replaced  by  water.  Conductivity  values  for  the 
bottom  sediment  and  frazil  deposits  are  the  same 
as  used  for  line  X6.  As  in  line  X6,  maximum  con¬ 
ductivity  values  occur  over  the  open  water  sections 
and  minima  occur  where  frazil  ice  exists  to  the 
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Figure  5.  Comparison  of  measured  EM31  data  for 
cross  section  X3A  (solid  line )  with  a  theoretical 
model  (dashed  tine). 


river  bottom.  The  modeled  and  observed  data 
agree  excellently  up  to  120  m,  especially  over  the 
water  channels  where  there  are  changes  in  at  due 
to  changes  in  both  solid  ice  and  water  depths.  Be¬ 
tween  120  and  180  m  the  agreement  is  mainly  qual¬ 
itative  as  there  are  apparently  more  changes  occur¬ 
ring  than  were  found  by  drilling.  Assuming  that 
the  value  1 .0  mS/m  is  still  correct  for  the  bottom 
sediments  in  this  section,  the  frazil  ice  conductiv¬ 
ity  would  have  to  be  raised  to  about  7  mS/m  to 
given  more  reasonable  quantitative  correlation  be¬ 
tween  measured  and  modeled  data  in  this  section. 
This  suggests  that  the  electrical  permeability  in 
this  section  of  frazil  ice  was  greater  than  in  other 
sections. 

X4 

The  observed  and  modeled  data  for  X4  are 
given  in  Figure  6.  There  is  good  reason  to  believe 
that  the  instrument  was  operating  incorrectly  for 
this  profile  as  there  is  less  quantitative  agreement 
between  observations  and  theory  than  before  over 
most  of  the  profile,  especially  over  the  main  chan¬ 
nel.  The  error  is  apparently  an  offset  of  about  1 
mS/m,  which  seems  less  between  60  and  95  m  be¬ 
cause  the  instrument  does  not  measure  below  0 
mS/m.  Nevertheless,  the  qualitative  agreement  is 
good  as  maxima  in  the  profile  occur  only  over  the 
open  water  sections.  The  maximum  seen  in  a,  in 
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Figure  6.  Comparison  of  measured  EM 3 1  data  for 
cross  section  X4  (solid  line )  with  a  theoretical 
model  (dashed  line). 
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the  dotted  profile  from  50  to  60  m  is  due  to  a  dra¬ 
matic  decrease  in  the  thickness  of  the  ice  cover  be¬ 
tween  the  40-  and  50-m  stations  (see  Appendix  B, 
Table  B3),  thus  allowing  the  more  conductive 
water  to  be  closer  to  the  instrument. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  conductivity  of  the  frazil  deposits  usually 
contrasted  well  with  the  conductivities  of  the  river 
water  and  the  solid  ice  cover,  bu‘  the  conductivity 
of  the  gravels  or  gravelly  sand  at  the  bottom  of  the 
Tanana  were  comparable  to  that  of  the  frazil  ice. 
Thus,  the  major  obstacle  to  magnetic  induction 
frazil  detection  would  be  shallow  streams  where 
bottom  sediments  could  be  confused  with  frazil 
ice.  Conductivity  values  far  below  those  expected 
for  solid  ice  over  water  would  be  a  significant  indi¬ 
cation  of  frazil  ice  deposits,  if  water  depth  were 
known  to  be  much  deeper  than  the  penetration 
depth  of  the  instrument.  Frazil  deposits  with  high 
volumetric  water  contents  will  be  more  conductive 
and  less  detectable.  Therefore  the  MI  method 
would  apply  primarily  to  dams  or  deposits  with 
low  water  content. 

Any  one  vertical  cut  in  a  cross  section  can  have 
as  many  as  four  different  layer  thicknesses  and 
five  different  conductivities.  The  thickness  of  the 
ice,  frazil  ice,  and  water  and  the  conductivities  of 
the  frazil  and  bottom  sediments  would  all  be  un¬ 
known.  Therefore,  the  use  of  the  MI  technique 
(without  drilling)  to  resolve  all  these  parameters 
requires  a  multimode  operation  in  which  either  the 
loop  orientation,  frequency,  and/or  coil  separa¬ 
tion  is  varied.  Use  of  one  instrument  in  the  vertical 
coplanar  (VCP)  mode  as  well  as  the  HCP  mode 
would  increase  sensitivity  to  frazil  deposits  but 
generally  decrease  sensitivity  to  depth.  The  use  of 
two  MI  instruments  with  different  coil  spacings 
would  supply  four  pieces  of  information,  enough 
to  solve  theoretically  for  the  truly  unknown  par¬ 
ameters  solid  ice  depth,  frazil  ice  or  water  depth, 
frazil  ice  conductivity,  and  bottom  sediment  con¬ 
ductivity,  assuming  either  water  or  frazil  ice  exists 
at  any  point  beneath  the  solid  ice  cover,  but  not 
both.  Thickness  of  the  ice  cover  can  be  measured 
separately  with  GPR  using  pulses  centered  above 
500  MHz,  which  has  demonstrated  sensitivity  to 
solid  ice  depth,  but  not  to  frazil  deposits  with  high 
water  content  (Arcone  and  Delaney,  in  prep). 
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APPENDIX  A:  DISCUSSION  OF  ERRORS 


The  sources  of  error  associated  with  the  MI 
method  include  internal  system  errors  and  external 
noise  errors.  System  errors  include  offset  error 
due  to  incorrect  calibration  and  zero  position  drift 
due  to  temperature  changes.  Noise  errors  mean  re¬ 
sponses  to  extraneous  electric  or  magnetic  Helds 
(coherent  or  random)  generated  by  artificial  or 
natural  sources.  Inaccurate  estimates  of  layer  con¬ 
ductivities  will  obviously  bias  any  modeling  at¬ 
tempt  at  interpretation. 

Offset  errors  are  nullified  by  calibrating  the 
zero  position  of  the  meter  over  an  area  of  known, 
low  conductivity,  preferably  less  than  10"4  S/m. 
This  procedure  was  repeated  several  times  during 
these  surveys.  Zero  position  drift  is  specified  by 
the  manufacturer  to  be  about  10  4  S/m  per  10°C 
change  in  temperature  for  the  quadrature  phase 
component.  This  is  a  minute  correction  and  the 
manufacturer  claims  reliable  operation  to  -40°C. 
All  surveys  were  performed  1  and  2  March  1986, 
during  which  time  temperatures  ranged  between 
-2  and  -35°F.  The  meter  contains  a  battery 
check  mode,  which  always  indicated  proper  func¬ 
tion.  Some  readings  made  on  land  during  this  time 
proved  inconsistent,  so  that  changing  water  or 
bottom  sediment  conductivity  are  not  believed  to 
be  a  possible  cause  for  the  difference  in  response 


levels  between  line  X4  and  the  other  two.  Prob¬ 
lems  related  to  temperature  are  therefore  still  sus¬ 
pect,  possibly  in  the  metering  system,  as  discussed 
below. 

Errors  due  to  magnetic  fields  generated  by  tel¬ 
luric  currents  are  possible  in  the  Arctic  due  to 
worldwide  thunderstorm  (sferics)  or  by  auroral  ac¬ 
tivity.  The  maximum  magnetic  Hux  density  gener¬ 
ated  by  the  instrument  transmitter  (a  coil  of  200 
turns  and  current  moment  of  about  1  amp-m1)  is 
about  4xl0‘7  Tesla,  a  value  easily  exceeded  by 
natural  fields.  However,  it  is  doubtful  that  suffi¬ 
cient  field  strength  exists  at  39  kHz  to  warrant  ser¬ 
ious  consideration  of  telluric  or  auroral  activity, 
which  is  generally  confined  to  frequencies  less 
than  10  Hz.  No  power  lines  or  underground  cables 
were  remotely  near  the  study  site.  The  most  plausi¬ 
ble  explanation  is  error  due  to  static  electric  fields 
generated  by  inadvertent  hand  contact  on  the 
meter  window.  This  problem  only  arises  in  very 
cold  temperatures  when  it  is  easily  induced  by  rub¬ 
bing  the  window.  Leakage  of  such  charge  appears 
to  be  slow.  The  best  solution  is  to  utilize  a  digital 
readout  system  (connected  to  a  data  logger), 
which  has  recently  become  available  from  the 
same  manufacturer. 
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APPENDIX  B:  MODELING  DATA 


The  following  tables  give  the  material  thicknesses  used  to  model  the  EM31  data.  An  addi¬ 
tional  (estimated)  91  cm  was  added  to  the  snow  thickness  to  account  for  the  operating  height 
of  the  instrument.  The  values  obtained  are  for  the  horizontal  coplanar  (HCP)  mode,  which  is 
also  referred  to  in  the  literature  as  the  vertical  magnetic  dipole  (VMD)  mode.  The  equation 
used  for  computing  apparent  conductivity  at  is  as  follows: 

a,  =  a, [1  —  /?(z,)]  -t-  <t2(/?(Zi)  —  /?(Zj)1  -t-  cr,[/?(Zi)  —  -t-  —  (Bl) 

where 


R(zJ  =  W(2zn/r)J  +  l  n  =  1,2,3....  (B2) 

The  quantities  an  are  the  successive  conductivities  of  the  layers,  zn  the  depths  from  the  sur¬ 
face  to  each  of  the  layer  interfaces,  and  rthe  intercoil  spacing,  which  is  3.66  m  for  the  EM31. 
This  layering  theory  may  be  applied  where  the  extent  of  the  layering  beneath  the  point  of 
measurement  extends  at  least  over  an  area  whose  radius  s  r. 


Table  Bl.  Snow,  ice,  and  water  layer  thicknesses  used  in  eq  Bl  to  derive  the  EM31  con¬ 
ductivity  values  (columns  6,  7)  for  line  X6. 


Dlst 

(m) 

Snow* 

(cm) 

Solid  Ice 
(cm) 

Frazil 

(cm) 

Water 

(cm) 

EM31 

(mS/m) 

Frazil 
replaced 
by  water 
(mS/m) 

10 

15 

30 

... 

_ 

0.80 

15 

15 

41 

— 

142 

3.22 

20 

15 

81 

— 

307 

4.22 

30 

15 

91 

— 

490 

4.89 

40 

15 

124 

— 

442 

4.20 

50 

15 

150 

353 

— 

0.47 

3.51 

60 

15 

124 

315 

— 

0.51 

3.65 

70 

15 

117 

231 

— 

0.54 

3.23 

80 

15 

135 

152 

— 

0.54 

2.43 

90 

15 

112 

61 

36 

0.98 

1.87 

100 

15 

112 

— 

163 

2.50 

110 

15 

97 

— 

213 

3.34 

120 

15 

71 

51 

91 

1.97 

2.92 

130 

15 

71 

— 

112 

2.55 

140 

15 

84 

— 

66 

1.81 

150 

15 

94 

— 

69 

1.79 

160 

15 

76 

— 

23 

1.14 

170 

15 

81 

8 

0.84 

*Not 

recorded . 

Values  In 

this  table 

are  the 

average  for  all 

other  lines. 

Table  B2.  Snow,  ice,  and  water  layer  thicknesses  used  in  eq  B1  to  derive  the  EM31  con¬ 
ductivity  values  (columns  6,  7)  for  line  X3A. 


Frazil 

replaced 


Disc 

(m) 

Snow 

(cm) 

Solid  ice 
(cm) 

Frazil 

(cm) 

Water 

(cm) 

EM31 

(mS/m) 

by  water 
(mS/m) 

15 

8 

71 

_ 

320 

4.56 

20 

5 

79 

— 

335 

4.57 

30 

2 

124 

— 

386 

3.95 

40 

1 

79 

— 

323 

4.37 

50 

2 

102 

— 

244 

3.51 

60 

1 

112 

— 

201 

3.09 

70 

10 

1 12 

170 

— 

0.57 

2.86 

80 

3 

99 

183 

— 

0.60 

3.19 

90 

10 

117 

155 

— 

0.57 

2.69 

100 

1 

122 

122 

— 

0.57 

2.27 

no 

23 

107 

155 

— 

0.57 

2.66 

120 

23 

79 

122 

43 

1.09 

3.01 

130 

30 

79 

71 

— 

0.63 

1.83 

140 

28 

69 

145 

— 

0.62 

2.85 

150 

25 

79 

119 

33 

0.99 

2.86 

160 

23 

91 

124 

— 

0.60 

2.46 

170 

30 

84 

99 

— 

0.61 

2.16 

180 

20 

71 

— 

81 

2.09 

190 

18 

91 

— 

33 

1.28 

200 

23 

64 

— 

48 

1.60 

210 

25 

76 

— 

1 

0.71 

Table  B3.  Snow,  ice,  and  water  layer  thicknesses  used  in  eq  B1  to  derive  the  EM31  con¬ 
ductivity  values  (columns  6,  7)  for  line  X4. 


Frazil 

replaced 


Dlst 

(m) 

Snow 

(cm) 

Solid  ice 
(cm) 

Frazil 

(cm) 

Water 

(cm) 

EM31 

(mS/m) 

by  water 
(mS/m) 

10 

3 

112 

_ 

348 

3.81 

15 

28 

74 

— 

442 

4.75 

20 

23 

74 

— 

686 

5.62 

30 

0 

84 

— 

922 

6.30 

40 

5 

145 

— 

688 

4.71 

50 

5 

107 

579 

— 

0.51 

5.08 

60 

5 

102 

561 

— 

0.51 

5.11 

70 

15 

94 

541 

— 

0.48 

4.60 

80 

20 

122 

406 

— 

0.49 

4.03 

90 

20 

99 

185 

— 

0.57 

3.02 

100 

23 

117 

— 

97 

1.95 

110 

23 

91 

— 

69 

1.76 

120 

23 

81 

— 

8 

0.83 

130 

23 

81 

— 

— 

0.69 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  repro¬ 
duced  below. 
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